Collapsing dynamics of attractive Bose-Einstein condensates 
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The self-similar collapse of 3D and quasi-2D atom condensates with negative scattering length is 
examined. 3D condensates are shown to blow up following the scenario of weak collapse: The in- 
ner core of the condensate diverges with an almost zero particle number, while its tail distribution 
spreads out to large distances with a constant density profile. For this case, the 3-body recombina- 
tion arrests the collapse, but it weakly dissipates the atoms. The confining trap then reforms the 
condensate at later times. In contrast, 2D condensates undergo a strong collapse: The atoms stay 
mainly located at center and recombination sequentially absorbs a significant amount of particles. 

PACS numbers : 03.75.Fi, 32.80.Pj, 42.65. Jx. 



A few years ago, Bose-Einstein condensates (BECs) 
were discovered in trapped clouds of alkali atoms ||],||. 
Among those, 7 Li atoms are known to be characterized by 
attractive interactions with a negative scattering length, 
so < |. This promotes a collapse- type instability, 
which yields a singular increase in the BEC wavefunc- 
tion. It occurs when the number of particles, N, exceeds 
a threshold value, N c ~ 'o/l a o|) where Iq = (h/muj) 1 ^ 2 
is the amplitude of zero point oscillations in the confin- 
ing trap of frequency to and m is the atom mass. BECs 
formed in 7 Li gas develop several sequences of collapse 
||. From the initial cooling, the condensate is first fed 
by the thermal cloud of uncondensed atoms. Then, once 
the number of atoms is above N c , the condensate sharply 
concentrates with an increasing density. As the density 
rises, inelastic collisions as 3-body molecular recombina- 
tion increase, which arrests the collapse. The thermal 
cloud is next believed to re-fill in the condensate, that 
reaches N > N c and then collapses again. The cycle 
of collapses repeats many times until the gas comes to 
equilibrium. This scenario is supported by numerical in- 
tegrations of the Gross-Pitaevskii (G.-P.) equation 

id t ijj = -A r tp + r 2 ip - \ip\ 2 ip - irj\ip\ m ^ + ijip, (1) 

where r = r/lo, t = tw/2, ip — -^(SttZq |a 1) 1 ^ 2 are 
the dimensionless time, coordinates and wavefunction 
of the BEC (tilde refers to physical quantities). Here, 
the loss/gain mechanisms are described by the last two 
terms of Eq. ([!]), where r\ <C 1 and 7 < 1 are the co- 
efficients related to the local recombinational decrease 
of the condensate density and to the flux of particles 
from the nonequilibrium thermal cloud to the conden- 
sate, respectively. At leading order, we consider 3-body 
recombinations with m = 4. The normalized number of 
particles, N — f \ip\ 2 dr, is related to its physical coun- 
terpart as N — iV(/o/87r|ao|) ~ 86. 7N. It was recently 
shown that stationary condensates defined by ij) s (r, t) — 
x(r, /u)e _vt with chemical potential /x are stable when- 
ever they satisfy dN eq /dfi < 0, where N eq (fi) = J \x\ 2 dr 



||. Along the curve N eq ([i), 3D condensates then re- 
quires N eq < N c = 14.45 for being stable. In physical 
units, this threshold yields the precise critical number 
for stability: iV c = 1252 atoms, for the 7 Li parameters 
ao = —1.45 nm and lo = 3.16 /im used in |0,||H]. 

Numerical simulations Q| of Eq. ([!]) revealed that, 
near the collapse moment, the condensate expands with 
a density profile exhibiting a low-amplitude, almost flat 
plateau, |-0| 2 r 2 — > const, from r > 1. This plateau-like 
behavior was again numerically found in where the 
free collapse was described as a "black hole". Follow- 
ing this scenario, the increase in the BEC density is fu- 
eled by particles drawn from throughout the outer region 
(r > 1), while the density outside form imploding ripples. 
However, it was suggested that the asymptotic plateau in 
r 2 \ip\ 2 may not be constant, but instead diverges in time. 
Thus, the singular dynamics of 3D collapsing BECs is 
still questionable and a self-consistent model for BEC 
collapse is actually missing. Describing the structure of 
such collapses is of utmost importance, in order to under- 
stand the influence of recombination, trap confinement 
and re-feeding over several collapse events. 

In this letter, we examine the self-similar nature of col- 
lapsing BECs. For 3D isotropic condensates the number 
of particles is analytically shown to vanish near center 
and, outside, a constant plateau in the density profile, 
r 2 |-0| 2 , is actually formed. From this dynamics, 3-body 
recombination arrests the collapse by removing a limited 
number of particles from the condensate. We also briefly 
investigate 2D condensates. As justified in [fr|, quasi- 
2D BECs can be produced from 3D atom clouds with 
a density frozen on a Gaussian shape oc e~ z ' °, when 
the particles are tightly confined along the longitudinal 
axis. In what follows, 2D and 3D condensates are con- 
sidered as isotropic, radially-symmetric objects, since the 
self-compression induced by the collapse dynamics make 
them have comparable sizes along each direction. 

Let us first discuss the inertial regime of collapse, for 
which we set 7 = r\ = in Eq. (nl). Collapse oc- 
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curs in the sense that the mean-square width of attrac- 
tive BECs, (r 2 ) = N- 1 J r 2 \ijj\ 2 dr, tends to zero m fi- 
nite time (see, e.g., Pitaevskii Q). From the inequality 
N < (r 2 ) f\Vip\ 2 dr @, the vanishing of (r 2 ) leads to 
the blow-up of the gradient norm, which in turn implies 
the divergence of the integral J \ip\ 4 dr in the conserved 
Hamiltonian of Eq. (m : 



H = J |V^| 2 (ir-i J |t/>| 4 dr + N(r 2 ). 



(2) 



By virtue of the mean- value theorem |?/;| 4 < 
maxrli/)! 2 J \ip\ dr, the maximum amplitude of the wave- 
function also blows up in finite time. It should be em- 
phasized that the blow-up generally occurs before (r 2 ) 
reaches zero [fl . To examine the shape of collapsing con- 
densates, we introduce the self-similarlike substitution: 



iP(r,t) = a- a (i)0(£T)e lAT - l/3 « 2/4 , 



(3) 



where £ = r/a(t), r(t) = L du/a 2 (u) and — — ah (dot 
means differentiation with respect to time). Here, the 
parameter A ~ —fj, must be positive in the absence of the 
trap for making <f> localized. The function a(t) represents 
the BEC scale length that vanishes as collapse develops, 
and 0(£, r) is a regular function with amplitude of order 
unity. As a(t) — > 0, it is assumed that <f> converges to an 
exactly self-similar form </>(£), which no longer depends 
explicitly on time, i.e., d T (j) — > 0. In this limit, the right 
balance between the two integrals in Eq. (Q) requires 
a = 1, in order to assure the finiteness of H. The particle 
number N = J \ip\ 2 dr then reads N = a D ^ 2 (t) f \(j)\ 2 d£ 
and the dynamics drastically changes following the space 
dimension number D. Setting thus a = 1 and plugging 
Eq. (0) into M) transforms the G.-P. equation into 



id T (f> + |0| 2 + e[^ 2 



ir/a 



(4) 



where A ? = C 1 " 13 ^^" 1 ^, e 
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i(/3 2 +/3 r ) and 



£^ = e _1 [A + i(3(D/2 — 1)] is viewed as a turning point. 
After an initial stage during which the trap gathers the 
particles at center, the wavefunction diverges freely on 
short time scales At < w -1 . It becomes hyper localized 
at r ~ with a(t) — > 0, so that the effect of trapping can 
be neglected. Near the collapse point, refeeding from 
the surrounding cloud is also inefficient. Moreover, 3- 
body recombination does not act, as long as the BEC 
radius satisfies rja 2 ~ m \<j>\ m <C |0| 2 . In this regime, Eq. 
(0) thus reduces to the self-similarly transformed nonlin- 
ear Schrodinger (NLS) equation, whose properties, accu- 
rately verified in P-pH, are recalled below. 

In the self-similar limit a(t) — > 0, d T <p — > 0, the time- 
dependent function e converges to /3 2 /4, which, for self- 
consistency, converges to a constant, eq. (ii) The so- 
lution (f> in Eq. (0) can be decomposed into a core 



contribution (f> c extending in the range £ < £t, i-e., 
r < tt = <z(t)|£r|, and a taiZ 4>t defined in the com- 
plementary spatial domain £ > £r, i-e., r > tt- (hi) 
The self-similar assumption <9 T </> — » holds provided that 
£ < £max = A\£T\/a(t), where A = const ^> 1. Self- 
similar solutions are thus limited by the boundary radius 
Tmax — A\£t\ ^> 1< Knowing this, the solution <j> de- 
composes as <f> = 4>c\t<£ T + 4>T\£ T <£<£ mam - In the long 
distance domain, the nonlinearity vanishes, so that </>t 
can be determined from the linear version of Eq. (Q) by 
means of WKB methods ||. As a result, the wavefunc- 
tion tp reads near the collapse point (e = /3 2 /4) as 



i>(r,t) = 



a(t) 



C{0) 



{M~,e)e-^^ \ O <_r<r T+{r/a)1+iX/0o 



\C(0)f 
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|r' T <r'<r max } 5 



(6) 
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f{f = - j + (- - l)(l + 21n2) - -arctan[^(- - 1)]. 

The scaling law a(t) in the inertial range of collapse is 
then identified through the continuity equation for i/>: 



d t \ip\ 2 r D ^dr = -2r I? -V| 2 5rarg(V')r 



-It) 



m+2 D-l 



r 
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dr + 27 



IVfr^-Vr, (7) 



where axg[ip(r = 0)] = 0. By applying the solution 
(H); (h) to Eq. (0) an( i Taylor-expanding C around 
e = /3 2 /4 = e 7 the contraction scale a(t) is indeed readily 
determined for 77 = 7 = from the dynamical system: 



£>-2, 



(8) 



where ci oc Re J 4>Qd e (f)\ eo dt; and C2 > are constants. 

We specify the structure of the collapse for 3D and 
quasi-2D BECs separately: 

1 - Three-dimensional condensates: Let us first de- 
scribe the inertial range of collapse. For D = 3, Eq. 
(Q) shows that rapidly attains a fixed point 0q > 
corresponding to the self-similar state T = 0. The 
scaling law a(t) ~ (t c — t) 1 / 2 follows, where t c de- 
notes the collapse moment. The characteristics of a 
3D collapse is that the number of particles is not pre- 
served self-similarly in the whole spatial domain, since 
N = a(t)N{<j>} with N{(/)} = J \(t)\ 2 dC By virtue of Eq. 
(p|), we see that, near the collapse instant, N behaves as 
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N ~ N core (t) + Ntau® = a(t)[JVT{0 c }(t) + iV{0r}(i)], 
where 

N{0 c } = 4n f T \<f> c \ 2 ed(i~0(l), (9) 
Jo 

A{<M = 4tt|C(/3)| 2 d£ ~ 47r|C(/3)| 2 r maa; /a(t), 

so that almost all A lies in the tail, N ta ii(t) — 
4:ir\C(f3o)\ 2 r max , as a(i) — * 0. Therefore, a 3D collapse 
takes place at the center of the trap where the wave- 
function = \cf>\/a(t) diverges. However, it is accom- 
panied by an expulsion of particles towards the large- 
distance domain r rr(i)- Thus, after the onset of 
collapse, the singularity develops, not by taking parti- 
cles from outside as in the "black hole" scenario pro- 
posed in ||, but by ejecting particles outward the core 
domain. This is a weak collapse, as originally defined 
in Ref. |Q for the free NLS equation. The solution ip 
blows up at center, near which the number of particles 
becomes zero, i.e., N core (t) = 4ir Jq T ^ \ip C ore\ 2 r 2 dr — > 0, 
as rr(t) ~ a(t) vanishes. Accordingly, the tail of the 
BEC wavefunction in the outer domain extends in space 
with the time-independent density: r 2 |?/>| 2 — > |C(/?o)| 2 
= const, deduced from Eq. (||). This constant de- 
pends on the values of /3q and \£t\- For the scaling law 
a(t) = aaitc — t) 1 / 2 with do set equal to the unity without 
loss of generality, /3q = 1/2 and numerical integrations of 
the 3D NLS equation with no trap (ll[| emphasize the 
values <fio = 0(0) = 1.39 and A = 0.545, which yield 
|£r| ~ 3.1 and |C(/3 )| 2 ^ 0.17 - 0.2. A tail contribu- 
tion r 2 \ip\ 2 /\ip(t = 0)| 2 about 0.2 in magnitude for r > 1 
seems compatible with the density profiles computed in 
Refs. in the vicinity of the collapse moment. The 
formation of stationary plateau-like density profiles was 
numerically confirmed in |l0|]ll| 1 by means of very fine 
numerical schemes that solved the rescaled equation (^) 
with high accuracy and could access huge growths in |-0| 2 . 
In contrast, the same plateau was claimed to diverge in 
time by only refining the spatial grid in ||. We suspect 
that numerical computations in Q suffered serious lack of 
resolution, which prevented the authors from concluding 
correctly on the constancy of r 2 |?/>| 2 at large distances. 

We now discuss the dissipative/gain regime, that in- 
volves collisional losses and feeding by the thermal 
cloud. The influence of these two effects on the con- 
densate particle number is described by the last two 
terms in the continuity equation (Q). By introduc- 
ing the self-similar shape (^), these terms are found to 
read ^rja 1 -" 1 J^ T |0 c | m+2 £ 2 d£ and 2-ya J^ ax |0r| 2 £ 2 d£ 
for a{t) <C 1, respectively. Hence, 3-body recombina- 
tion mainly acts on the core part of the solution \%j) c \ = 
\(j) c \/a(t), whereas re-feeding by the surrounding cloud is 
efficient in the outer region, where the main amount of 
particles is residing, with 2j J rmax \ip\ 2 r 2 dr ~ 2jN ta ii{t). 



Collisions begin to be active when a(t) decreases so 
much that rja 2 - m \(j)\ m ~ \(j)\ 2 [Eq. (|)], i.e., this contri- 
bution saturates the blow-up induced by the cubic non- 
linearity, for m > 2. By means of this relation, the num- 
ber of particles ANi oss lost from the condensate during 
one collapse event is then yielded by that in the core re- 
gion, N core . Inelastic collisions are thus able to stop the 
collapse, but they cannot remove a lot of particles, the 
major part of A ~ N ta ii starting to be transferred to 
large distances via the weak collapse. Explicitly, we find 
by using Eq. (^) expressed with arbitrary m and D: 

AN loss ~ -2 D -\M£,T\) D Kri {D - 2)l(m - 2) lD. (10) 

For D = 3 and m = 4, the particle loss is estimated 
as ANi oss — ^y/v^Mr/^ — 0.5, when we choose the 
value 77 = 2.2 x 10~ 6 physically justified in |J and em- 
ploy 0o = 1-39 and £t = 3.1. Thus, in physical units, the 
number of atoms removed during a single blow-up event 
is ANioss = 44 atoms. This loss of particles agrees with 
recent numerical observations [fl3|| . Here, approximately 
300 particles were lost from a condensate of initial num- 
ber Ao = 1260 atoms, in a sequence of 5-6 individual 
blow-up events with peak amplitudes IV'Imaa: ~ 7 y 1 / 2 -m 
As discovered by Vlasov et al. and confirmed in [jll| , 
an untrapped, 3D collapsing field is damped within a 
very short time window, in which the field experiences 
secondary blow-up events of duration ~ ^ 2 / m ~ 2 ; as long 
as nonlinear dissipation remains active. These arise while 
most of the atoms form a broad disribution surrounding 
the collapsing core and are directed outwards. Collapse 
is not arrested abruptly, but it is distributed through a 
series of individual spikes. This "distributed collapse" 
p4| resembles the "intermittent implosions" numerically 
detected in p3| . Note that these intermittent implosions, 
which overlap as one collapse sequence over a trap period, 
may not be resolved in real experiments. Their integrated 
effect should only yield a smooth decrease in A. 

Once collapse is limited, it is necessary to under- 
stand the mechanism underlying reformation and peri- 
odic resurgences of spiky amplitudes in the condensate 
||. We emphasize that after the collapse stage most 
particles lie in the outer domain 1 < r < r max , where 
the atoms feel both the trap curvature and re-feeding 
by the thermal cloud. The latter re-injects particles 
into the BEC according to N(t) ~ A ta i;e 27 *, i.e., on 
time scales t ~ (2^f)^ 1 In (1 + AA/A ). For reasonable 
values of 7 = 3x 10 ~ 3 ||, the physical time for re- 
filling the condensate with AN/Nq ps 0.25 is too long 
for justifying the spiky oscillations occurring periodi- 
cally at times ~ w^ 1 , mentioned in, e.g., |l3). We have 
thus to analyze the stage of trap reconfinement at ear- 
lier instants for which refeeding is negligible. To this 
aim, we define the centroid of the expanded structure by 
X(t) = A -1 j r\ij]\ 2 dr. From algebraic manipulations of 
Eq. (|I]) with 7 = 77 = 0, this centroid is found to obey 
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the relation X + AX = 0. Simple phase transformation 
|l5f keeps equation (|l|) unchanged when passing over to 
the frame moving with X(t). It is then easy to conclude 
that, from a reference coordinate X(to) ^ that locates 
the BEC expanded in the outer domain at a given time 
to, the condensate will return to the center of the trap 



over physical times = 7r/2w 



after the collapse mo- 



ment. Once reconfmed with N < N c , the condensate is 
close to an equilibrium, that is described by the ground 
state solution ip s = x(r, /i)e _v * of the G.-P. equation. 
As shown in Ref. 0] (see also ground states are 

stable with dN eq /d/i < for rather broad condensates 
(0.72 < n < 3). They are unstable and collapse with 
dN eq /d[i > for rather narrow condensates (/i < 0.72). 
In this range, the threshold N c > N eq decreases with 
as N c {ijl) < 18.94/|/x|. It is thus highy problablc that un- 
der the strong compression induced by the trap, the equi- 
librium state confined at center corresponds to a narrow, 
unstable ip s undergoing a new collapse sequence with a 
lower critical number of particles. This may repeat over 
several trap periods, until the refeeding becomes a key- 
player for making the resulting BEC bifurcate to a stable 
equilibrium reached at a weak number of atoms [0] . 

2- Quasi-two-dimensional condensates: The dynamics 
significantly changes when one considers quasi-2D BECs. 
For r\ = 7 = 0, the scaling law a(t) inferred from Eq. 
(0) with D = 2 indeed behaves with a twice-logarithmic 
correction: a{t) ~ a Q y/t^t/ ^lnln [\/{t c - t)] 
Since (3 ~ 7r/lnr(t) -> as r(t) ~ ln[l/< c -t] -> +oo, 
the exponential contribution in the tail amplitude of Eq. 
(0) decreases to zero, while the core boundary £ T = 
increases slowly to infinity when t — > t c . As a result, 
collapse takes place with a core solution <f> c providing 
the principal contribution in the wavefunction tp. N is 
mainly given by J \<j> c \ 2 d£,, which relaxes to the critical 
value N c = j \x\ 2 d£, — 11.7. Below this number, sta- 
tionary condensates are stable Q. So, most of parti- 
cles stay located around the center of the trap, when 
collapse develops. This scenario meets the definition of 
the so-called strong collapse, which contrasts with the 
3D weak collapse that promotes a leak of atoms to large 
distances. The relative variations in the particle num- 
ber due to collisional losses are governed by the estimate 
AN loss ~ -2</>Q V 2 ~ m / \(t)c\ m+2 d-C As a(t) decreases, 
the part burnt by recombination is here estimated by 
ANi oss /N ~ —2/(1 + m/2), when we apply the evalu- 
ation na 2 ~ m « |</>| 2_m and model 4> around a Gaussian 
with intensity Iq: 4> — y/loe~^ I" 1 . For to = 4, 3-body 
recombination then removes about more than one-half 
of the condensed atoms. Although rough, this estimate 
does not explicitly vary with the coefficient r\, as already 
indicated by Eq. (10), and it slightly decreases with to, 
in agreement with the numerical observations of Ref. [[T"l| . 
Collapse is thus arrested by 3-body recombinations that 
sequentially burn a substantial amount of atoms per col- 



lapse event. This agrees with the numerical results of 
Refs. from which AN loss is about 0.25-0.3 x N c . 

In conclusion, we have shown that BECs with attrac- 
tive interactions collapse with a mean radius contract- 
ing like y/t c — t at leading order, where t c denotes the 
collapse moment. In 3D, the collapse is weak: the am- 
plitude of the wavefunction blows up at center, while 
particles are expelled outwards with a constant density 
profile r 2 \ip\ 2 — > const. Recombination losses limit the 
collapse, but they only damp a few percents of the atoms 
at each blow-up. Several blow-up events develop within 
one collapse sequence. The condensate is then reconfmed 
by the trap and can further undergo more collapse cycles, 
as long as a stable state has not been reached. Quasi-2D 
condensates are subject to a strong collapse, in the sense 
that the number of particles remains mostly confined at 
center. In this case, recombination removes a significant 
fraction (up to 0.5) of particles per collapse event. 
This work was partly supported by the Danish Natural 
Sciences Foundation ( |snf-grant 9903273 ). 
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